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Abstract 

We are interested in modeling Darwinian evolution resulting from the in- 
terplay of phenotypic variation and natural selection through ecological in- 
teractions. The population is modeled as a stochastic point process whose 
generator captures the probabilistic dynamics over continuous time of birth, 
mutation, and death, as influenced by each individual's trait values, and in- 
teractions between individuals. An offspring usually inherits the trait values 
of her progenitor, except when a random mutation causes the offspring to take 
an instantaneous mutation step at birth to new trait values. In the case we are 
interested in, the probability distribution of mutations has a heavy tail and 
belongs to the domain of attraction of a stable law. We investigate the large- 
population limit with allometric demographies: larger populations made up 
of smaller individuals which reproduce and die faster, as is typical for micro- 
organisms. We show that depending on the allometry coefficient the limit 
behavior of the population process can be approximated by nonlinear Levy 
flights of different nature: either deterministic, in the form of nonlocal frac- 
tional reaction-diffusion equations, or stochastic, as nonlinear super-processes 
with the underlying reaction and a fractional diffusion operator. These ap- 
proximation results demonstrate the existence of such nontrivial fractional 
objects; their uniqueness is also proved. 
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1 Introduction 



In this paper, we are interested in modeling the dynamics of populations as driven 
by the interplay of phenotypic variation and natural selection operating through 
ecological interactions, i.e.. Darwinian evolution. The population is modeled as 
a stochastic Markov point process whose generator captures the probabilistic dy- 
namics over continuous time of birth, mutation and death, as influenced by each 
individual's trait values and interactions between individuals. The adaptive nature 
of a trait implies that an offspring usually inherits the trait values of her progenitor, 
except when a mutation occurs. In this case, the offspring makes an instantaneous 
mutation step at birth to new trait values. This microscopic point of view has been 
heuristically introduced in Bolker-Pacala [5] and Dieckmann-Law [12]. It has been 
rigorously developed first in Fournier-Meleard [15] for spatial seed models and by 
Champagnat-Ferriere-Meleard [8J, [9J for phenotypic trait structured populations 
when the mutation kernel behaves essentially as a Gaussian law. The aim in this 
work is to study the case where a mutant individual can be significantely different 
from his ancestor. More precisely, the mutation kernel will be assumed to have a 
heavy tail and to belong to the domain of attraction of a stable law. If the traits de- 
scribe a spatial dispersion, for instance for seeds, we thus assume that seed offsprings 
instantaneously jump far from the mother seed because of availability of resources or 
wind. In the case of phenotypic traits, the heavy tail mutation assumption says that 
if a mutant offspring is too close to the mother trait, then it is so deleterious that it 
cannot be observed. A trait which quantifies the aggressivity level is an example of 
such a situation. 

In the context of ecology several authors have considered the fractional reaction- 
diffusion model and we would like to acknowledge some of their contributions here. 
In particular, Baeumer, Kovacs, and Meerschaert [2] considered fractional repro- 
duction-dispersal equations and heavy tail dispersal kernels. Their paper also con- 
tains an exhaustive review of the literature on the population spreads with extreme 
patterns of dispersal and reproduction. Another population djTiamics model in- 
troduced by Gurney and Nisbet [TB] led to strongly nonlinear partial differential 
equation of the porous media type; we have studied a probabilistic fractional frame- 
work for related stochastic and partial differential equations in [19] . 

Each individual is characterized by a real- valued trait x describing a phenotypic, 
or a spatial parameter. The birth and death rates of this individual depend on its 
trait X and also on the environment through the influence of the others individuals 
alive. In case the offspring produced by the individual with trait x carries a mutated 
trait, then the difference between this mutated trait and x is distributed according 
to M{x, dh) such that 
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with a G (0, 2) not depending on x. The inclusion of the dependence of the mutation 
rate on the trait itself, as expressed in the dependence of M, and C, on x, is an 
essential part of our model even though it considerably complicates the mathematical 
framework of this paper. 

An example for M{x, dh) could be the Pareto law with density f l{|/i|>i}/|^|"'^^", but 
a mutation law equal to the Pareto law (up to a constant) outside a given interval 
of the form [—a, +a], and constant inside, is also covered by our model. The latter 
example corresponds to a distribution of mutant traits which is uniform in a small 
neighborhood around the mother's trait, but decreases, with heavy tail for more 
distant traits. 

We investigate the large population limits, K ^ oo, with allometric demographies: 
larger populations made up of smaller individuals which reproduce and die faster. 
This leads to systems in which organisms have short lives and reproduce fast while 
their colonies or populations grow or decline on a slow time scale. Typically, these 
assumptions are appropriate for microorganisms such as bacteria or plankton. The 
allometric effect will be modeled by a dominant birth and death rate of order K^, 
T] > 0. 

In [H] , such asymptotic approximations have been studied in the case of small muta- 
tions with a Gaussian, thin tail distribution. In such a situation, for the allometric 
exponent r] < 1, and large enough population, the individual population stochas- 
tic process is approximated by the solution of a classical deterministic nonlinear 
reaction-diffusion equation. In our present heavy tail setting we will prove that in a 
similar limit, the population process is approximated by the solution of a nonlocal 
nonlinear partial differential equation driven by a fractional Laplacian and involv- 
ing a reaction term. Note that, as a byproduct, this result provides a proof of the 
existence of weak solutions of such nonclassical equations. Separately, employing 
purely analytical techniques, we will also prove the uniqueness of their solutions. 
The case of the allometric exponent rj = 1 is significantly different as stochasticity 
remains present in the limit due to the demographic acceleration. The population 
process is then approximated by a nonlinear super-process with underlying reaction 
and fractional diffusion. 

Our work provides a rigorous derivation of macroscopic models from microscopic 
dynamics (hydrodynamic limit) for a large class of nonlinear equations involving 
Levy flights which naturally appear in evolutionary ecology and population dynam- 
ics. Our use of probabilistic tools of interacting particle systems is essential and 
provides the information about the scales (between mutation amplitude and popu- 
lation size) at which such models are justified. 

It should be mentioned at this point that fractional reaction-diffusion equations have 
been suggested and studied as models for several physical phenomena. Thus, Del- 
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Castillo-Negrete, Carreras, and Lynch [TT] investigated front dynamics in reaction- 
diffusion systems with Levy flights, and Henry, Langlands and Wearne [17] , studied 
Turing pattern formation in fractional activator-inhibitor systems. Also, Saxena, 
Mathai, and Haubold found explicit solutions of the fractional react ion- diffusions 
equations in terms of the Mittag-Leffier functions which are suitable for numerical 
computations; they also considered the situation where the time derivative is also 
replaced by a fractional derivative of order less than one. Another numerical method 
for finding solutions of such equations can be found in Baeumer, Kovacs, and Meer- 
schaert [3J. 

This paper starts (Section 2) with the description of the reproduction and death 
mechanisms for the individuals of the population we are interested in. The main 
convergence results based on a large population limit are stated. Thus Theorem 13.21 
shows that an allometric effect of order K^, with rj e (0, 1), leads to a determinis- 
tic, nonlinear integro-differential equation driven by a nonlocal fractional Laplacian 
operator independent of t], while in Theorem 13.31 we show that an allometric effect 
of order K {rj = 1) yields, in the large population limit, a stochastic measure- valued 
process depending on the acceleration rate of the birth-and-death process. Thus 
(demographic) stochasticity appears as the allometric exponent takes on the value 
1] = 1. Also, Prop osit ion 13 . 1 1 est ablishes that the mutation kernel conveniently renor- 
malized behaves approximately as a jump kernel with the heavy tail jump measure 
a{x)dz/\z\^~^°' . Section 4 contains technical lemmas needed in the proof of the main 
results of Section 3. In particular we clarify a key technical point that was eluded in 
[T5] . [H], allowing to deduce the tightness of the measure valued population process 
for the weak topology from the tightness for the vague topology (see Lemma and 
Remark 14.41 and Step 2 of the proof of Theorem 13. 21) . Section 5 contains proofs of the 
main theorems stated in Section 3 using the measure-valued martingale properties 
of the population process. 



2 Population point process 

As in [9], the evolving population is modelled by a stochastic system of interacting 
individuals, where each individual is characterized by a phenotypic trait. This trait 
is described quantitatively by a real number. We assume that the parameter K 
scales the initial number of individuals. To observe a nontrivial limit behavior of 
the system as K grows to infinity it is necessary to attach to each individual the 
weight Thus our system evolves in the subset A^x of the set Mp of finite non- 
negative measures on M consisting of all finite point measures with weight 

Here and below, 5^. denotes the Dirac mass at x. For any m G Mp, any measurable 
function / on M, we set (m, /) = J^fdm. Following [S], we describe the population 
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by 

JK 

-" = :^E'^-.- (2-1) 

i=l 

with I I" G N denoting the number of individuals ahve at time t, and X^, 
describing the individuals' traits (in M). 

The population measure- valued process evolves as a birth and death process with 
mutation and selection. More precisely, an individual can give birth or die. The 
death can be natural, or can be due to the competition pressure exerted by other 
individuals (for instance, by sharing food). At birth, the offspring can inherit the 
trait of its parent, or can mutate to another trait with some positive probability. 
We assume that the population is of order K and that for each individual, birth 
and death occur at the rate of order K^, for some < < 1, while preserving the 
demographic balance. More precisely, the main assumptions on the birth and death 
dynamics are summarized below. 



Scaling Assumptions: For a population = Xl[=i ^x^j ^-nd a trait a; G M, 
the birth and death rates are scaled with the system's size according to the 
following rules: 



bK{x,u) = K\{x) + b{x,^^V{x-x'j)=K"r{x) + b{x,V*u{x)), 

1=1 

dK{x,u) = ii'''r(a;) + (i(a;, — ^^/(x-a;*)) = /r''r(a;) + (i(a;,t/*t/(a;)), 

1=1 

(2.2) 

where b and d are continuous functions on M^, and * denotes the convolution 
operation. The allometric effect (smaller individuals reproduce and die faster) 
is parametrized by the exponent r) and a trait-dependent function r which is 
assumed to be positive and bounded on M. 

Assumptions (HI): The interaction kernels V, and U, affecting, respec- 
tively, reproduction and mortality rates, are continuous functions on R. In 
addition, there exist constants f, b, d, U, and V, such that, for x,z eM, 

< r{x) <f ; <b{x,z) <b, (2.3) 
< d{x,z) <d{l + \z\), (2.4) 
< U(x) <U ; 0< V{x) < V. (2.5) 

Note that the death rate d is not assumed to be bounded but its growth in 
variable z is at most linear. This dependence models a possible competition 
between individuals (e.g., for shared resources), increasing their death rate. 
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Assumptions (HI) ensure that there exists a constant (7 > 0, such that the 
total event rate for a population counting measure — ^ ^i=i obtained as 
the sum of all event rates, is bounded by (7/(1 + 7) , where I is the population 
size. 

Assumption (H2): When an individual with trait x gives birth, it can pro- 
duce a mutant offspring with probability p{x). Otherwise, with probability 
1 — p{x)^ the offspring carries the same trait x as its ancestor. If a mutation 
occurs, the mutated offspring instantly acquires a new trait x + h, where h 
is selected randomly according to the mutation step measure Mk{x, dh). We 
only consider mutations which have heavy tail distributions. More precisely, 
we assume that the probability measure Mk{x, dh) is the law of the random 
variable 

X(x) 



and a bounded cr : IR — >■ 1R+, 



where X{x) is a symmetric random variable such that, for some a e (0, 2), 

0. (2.6) 



lim sup 



a 



As an immediate corollary of Assumption (H2) we have 

2a{x) 



/ MK{x,dh) 

J\h\>v 



'\h\>y yK <.^oo a K^ya^ 

and, for a e (1, 2), also 

hMK{x,dh) = 0. 



Example 2.1. One can choose the random variable X[x) with Pareto's law inde- 
pendent of a;, that is, with the density 

a 1 



Example 2.2. Another possibility is to take X{x) with the Parcto law (up to 
multiplicative a constant) outside a given interval of the form [—a, +a], and constant 
inside. This choice corresponds to a distribution of mutant traits which is uniform 
in a small neighborhood around the mother's trait, but decreases, with heavy tail, 
for more distant traits. 
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We refer to Fournier-Meleard [15] or Champagnat-Ferriere-Meleard [8] for a pathwise 
construction of a point measure-valued Markov process {h'f)t>o satisfying Assump- 
tions (HI). The infinitesimal generator of its Markovian dynamics is given, for each 
finite point measure z/, by the expression 

L^Uiu) = K / (1 - p{x)) (K^rix) + b{x, V * z/(x)))(0(z/ + ^6^) - <p{v))v{dx) 
Jr K 

+ K [ p{x) {K'>r{x) + b{x,V*iy{x))) [ {<j){iy + ^6^^^) - <f>{'y))MK{x,dh)iy{dx) 
Jr Jr 

+ K [ (K^rix) + d{x, U * z/(x)))(0(z/ - —6^) - (p{u))u{dx). (2.7) 
Jr K 

The first term of fl2.7p captures the effect of births without mutation, the second 
term that of births with mutation, and the last term that of deaths. 

If E{{i'q , 1)^) < +00, then for any T < oo, £'(supjg[o,T] (^/^? 0^) ^ ^ (^^^ Lemma 
14. ip . Thus, for any measurable functions on Mp such that \ + |L^0(z/)| < 
C*(l + (z/, 1)^), the process 

<P{u^) - <l>{4) - fL^<P{uf)ds (2.8) 
Jo 

is a martingale. In particular, in view of f l2.7p . for each measurable bounded function 
/, 

= (z/f , /) - (4, f) 

f*t /• 

(6(x, V * vf{x)) - d(x, U * uf{x)))f{x)uf{dx)ds (2.9) 

^ p(x) (KV(x) + 6(x, F * z/f (a;))) f / /(x + h)MK{x, dh) - f{x) ) vf{dx)ds, 
is a square integrable martingale with quadratic variation 

(M^'/), = i-/ /" ! {2K'^r{x) + b{x,V*uf{x))+d{x,U*uf{x)))f{x)uf{dx)ds 
K [JoJr 

t r 



+ / p(a;) (irV(x)+6(x,y*z/f (x))) / f{x+h)MKix,dh)-f\x))iy^{dx)ds 
Jo Jr \Jr ) . 

(2.10) 

3 Scaling limits of population point processes 

Our aim is now to make the population size's scaling K tend to infinity, in a scale 
accelerating the births and deaths, and making the mutation steps smaller and 
smaller. We begin by explaining why a fractional Laplacian term appears in the 
limit, K — )■ oo, of the last term of the r.h.s. of (12. 9p . 
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Proposition 3.1 Under (112) , 

(i) ifae (1, 2), then for any f e Ci{R), J^{f{x + h) - f{x))MK{x, dh) 
converges to a{x) J^{f{x + z) — f{x) — f'{x)zl{\z\<i}) \z\^+c, uniformly for x e R as 
K oo. 

(a) if a e (0, 1], then the same result holds for f e C'^(R) compactly supported. 

Proof . By Fubini's theorem and since f{x + h) — f{x) = f'{x + z)dz, 

f l{h>o}ifi^ + f{x))MK{x, dh) = f f\x + z)¥{X{x) > Kiz)dz. 
Jr Jo 



Treating in the same way the integral for /i < and using the symmetry of X{x), 
one deduces that 

^'^'^'^-J'^'-'^ n\X{x)\ > Klz)dz. 



I {f{x + h)-f{x))MK{x,dh)^ f 
Jr Jo 



By integration by parts, 

dz /■+°° dz 

{fix + z)- fix -z))- = a {fix + z) + fix -z)- 2fix)) 



l+a 

dz 



= a / ifix + z) - fix) - f'ix)zl{\^\<i^)- 

Jr 

Combining both equahties one deduces that 

R(x)'^K'' j (f{x + h)-f{x))MK(x,dh)-a(x) j (/(i + 2) - /(x) - /'Wzl,|.|<„) 
- /'(. + .) - fix - z) ^^,,„p(|^(^)|"> ^J,) _ 2.(x)^ 



a I z 



Therefore 

< liriloo sup «"P(|X(X)|>«)-^ 
u>G,xeR Oi 



+ ||r||oo sup M"P(|X(a;)| >m)- 



u>K^,xeR 



2aix) 



+ sup 

u>Ka,xeR 



u"Fi\Xix)\ > u) 



2aix) 



a 



a 

»+oo 



z^-'^dz 



z^-'^dz 



2ct{x) 



By (H2) sup^^^^_^^^ u-n\Xix)\ > u) - ^ 
for V > 0, 

2aix) 



/■+°° dz 
I \f'ix + z)-f'ix-z)\ — . 

(3.1) 

tends to as — )■ 00. Moreover 



sup 

u>0,x&R 



u^'FilXix)] > u)- 



a 



2 2aix'\ 
< 't;°+- super (x)+ sup m"P(|X(x)| > ^ 



with the right-hand-side finite when v is large enough. 

Last, /+°° If'ix + z)- fix -z)\^is smaller than ||/'|U f if a > 1. If 
a G (0, 1], it is smaller than ||/'||oo Ji~^^^' fl when the compactly supported function 
f{y) vanishes for \y\ > C. Hence (13.11) implies the desired uniform convergence as 
K oo. 

I 

In what follows, we will denote by a the product function defined by 

a{x)=p{x)r{x)a{x) 

and by a the function 

a{x) = a^^°{x) = {p{x)r{x)a{x))^^" . 

The nature of the scaling limit of the population point processes strongly 
depends on the value of the allometric exponent rj. For < r/ < 1, the limit is 
deterministic and our convergence result is described in the following theorem: 

Theorem 3.2 (i) Suppose that Assumptions (HI) and (H2) are satisfied, < r] < 1 
and the product function a is continuous. Additionally, assume that, as K oo, 
the initial conditions Uq converge in law, and for the weak topology on Mp, to a 
finite deterministic measure one? that 

supE((z/o^,l)^) < +00. (3.2) 

K 

Then, for each T > 0, the laws (Q^) of the processes (z/^) in B){[0,T], Mp) (with 
Mp endowed with the weak convergence topology) are tight. Moreover, the weak limit 
of each of their convergent subsequences gives full weight to the process {^t)t>o € 
C{[0,T], Mp) satisfying the following condition: for each function f G C^(M), 

(6, /) = (eo, /) + / / {bix, V * U^)) - dix, U * U^)))f{x)Udx)ds 
Jo Jr 

^ fo L ^^""^ ( L^^^"" ^^""^ ~ /'(^)^l{|'^l<i})^) Udx)d.s. (3.3) 

(ii) If we assume additionally that b{x,z), and d{x,z), are Lipschitz continuous 
in z, uniformly for x G M, and that a is Lipschitz continuous, then (13.31) has at most 
one solution such that sup^gjQ {C,t, 1) < +oo, and, as K ^ oo, the processes (z/^) 
converge weakly to this unique solution. 

(Hi) Finally, if cr E Cf, i.e., it is bounded together with its derivatives of order 
< 3, and the product function a{x) > 0, for all x G M, then, for each t > 0, the 
measure has a density with respect to the Lebesgue measure. 
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The following remarks are immediate corollaries of the above Theorem. 

Remark 3.1. The limit (13.31) does not depend on 77 G (0, 1). As will appear in the 
proof, this is implied by the fact that the growth rate — does not depend on 
7], and that the mutation kernel Mk{x, z) compensates exactly the dispersion in the 
trait space induced by the acceleration of the births with mutations. 

Remark 3.2. In the case considered in Theorem 13.21 (in), Eq. f l3.3p may be written 
in the form 

dtUx) = [h{x, V * Ux)) - d{x, U * + D''{a^t){x), (3.4) 

where 

DVi^) = [ (fix + h)- fix) - /(x)M||,|<i|)-^, 
JR 

denotes the fractional Laplacian of order a. Note that, by the change of variable 
z = h/a{x), 

a(x)D"/(x) = / (fix + aix)z) - fix) - f(x)a(x)^l||,|<i|)— — . (3.5) 
Jr 1^1 



Remark 3.3. Theorem 13.21 also proves the existence of a weak solution for (13. 4p . 
Equation (13. 4p generalizes the Fisher reaction-diffusion equation known from classi- 
cal population genetics (see e.g. |6j) with a fractional Laplacian term replacing the 
classical Laplacian. It justifies the Levy flight modeling in ecology as approximations 
for models in which mutations have heavy tail. Such models were judged justified by 
real-life data in many areas of physics and economics, and in some popular literature 
devoted to the 2008 global financial crisis. 

In the case of the allometric exponent rj = 1 the scaling limit of the population 
point processes z/^ has a richer structure of a nonlinear stochastic superprocess 
which is described below. 



Theorem 3.3 (i) Suppose that Assumptions (HI) and (H2) are satisfied, 77 = 1 and 
that a is continuous. Additionally, assume that, as K 00, the initial conditions 
Vq converge in law, and for the weak topology on Mp, to a finite (possibly random) 
measure Xq, and that 

supE((z/o^,l)^) < +00. (3.6) 

K 

Then, for each T > 0, the laws of the processes G D([0, T],Mj7) are tight and 
the limiting values are superprocesses X G C([0, T], Mp) satisfying the following two 
conditions: 

sup E{{Xt,l)^) <oo, (3.7) 
te[o,r] 
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and, for any f G (M), 

Mi = {Xt, f) - (Xo, f) - [ [ H^) Dy{x)X,idx)ds 

Jo Jm. 

-II f{x){b{x,V*X,{x))-d{x,U*X,{x)))Xs{dx)ds (3.8) 
Jo Jr 

is a continuous martingale with the quadratic variation 

{Mf)t = 2 I I r{x)f{x)Xs{dx)ds. (3.9) 
Jo Jr 

(a) Assume moreover that a is Lipschitz continuous and that r is bounded from 
below by a positive constant. Then there is a unique such limiting superprocess. 

Remark 3.4- Here again, as in the case of Theorem 13.21 Theorem 13.31 yields the 
existence of a measure-valued process X which is a weak solution of the stochastic 
partial differential equation with fractional diffusion operator 

dtXtix) = (^bix, V * Xtix)) - dix, U * Xt{x))^Xtix) + D''{dXt){x) + M^, 

where Mt is a random fluctuation term reflecting the demographic stochasticity of 
this fast birth-and-death process; the process is faster than the accelerated birth-and- 
death process which led to the deterministic reaction-diffusion approximation (13.41) . 

4 Auxiliary Lemmas 

In this section we provide auxiliary lemmas needed in the proofs of Theorems 13.21 
and 13.31 The latter will be completed in Section 5. We begin with a lemma which 
gives uniform estimates for the moments of the population process. Its proof can be 
easily adapted from the the proof of an analogous result in p!5] and is thus omitted. 

Lemma 4.1 Assume that p > 2, and supj^ i?((i/(f" , 1)^) < +oo. Then 




To provide mass control for z/ , we first need to study the action of the fractional 
Laplacian on functions 

/„(a;) =V^(OV(|x|-(n-l)) Al), n G N, 

where the function ip{x) = 6x^ — 15x^ + lOx'^ G C^. Observe that ip is non- decreasing 
on [0, 1], and such that 

ijj{0) = ip'iO) = V"(0) = 1 - tpil) = = ^"(1) = 0. 
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Lemma 4.2 For each n E N* , the function fn is in C"^ , even, non- decreasing on 
M+, equal to on [—{n — — 1], and to 1 on {~n,nY. In particular /q = 1. 
Moreover, 

sup \D"fn{x)\ < +00, 

neN*,a;eR 

and, for each n eW, and x G {—n + 1, n — 1), 

2 

a{n — 1 — |a;|)" 
Finally, under Assumption (H2), 



\D''fn{x)\ < 



lim / + h)- fn{x))MK{x, dh) = a{x)D''fn{x), 



uniformly, for n eN* , and x G 



Proof . For n E W, let us check the three last statements for /„, the others 
being obvious. Using the Taylor expansion for |?/| < 1 and remarking that, for each 
x,z eR, \fn{z) - fn{x)\ < 1, oue has 



\D''fnix)\ < 



sup,emWiz)\y^ dy 



X 



|2/I<1 



\y 



l+a 



+ 



dy 



\y\>i \y 



l+a 



< 



sup^6[o,i] 2 



2 — a a 

For X E (— 72 + 1, n — 1), since = /^(x) = 0, and fn{x+y) = if |?/| < n — l — \x\ 



l^"/n(x) 



fn{x + y)dy 



I l+a 



< 



dy 



H>n-i-|x| a(n-l-|x|)"' 



Since sup„gp^* .^^^ l/n(a;)l = sup^.g[o i] Wix)] < +00, writing for the function 
/„ we see that to establish the uniform convergence we only need to check that 
sup„eN*,xeM/i^°° Ifni^ + ^) - f'ni^ " < Since sup^gp^.^^gR | = 

sup^.g[o,i] < +00, and vanishes outside [—n, —n + 1] U [n — 1,^], 



sup 

neN*,xeK Jl 



dz 

\fn{^ + z)-fA^-z)\—< sup W{X)\ 



2^° 



xe[o,i] 



dz 



which concludes the proof. 

The next lemma provides mass control for the sequence {v^)- 
Lemma 4.3 Under the assumptions of Theorem \3.2\ or Theorem \3.3[. 



lim lim sup ( sup(z/j , /„) = 0. 
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Proof . Let M/^'" denote the square integrable martingale defined by (12.91) with /„ 
replacing /. The boundedness of r and a together with the last assertion in Lemma 
14.21 ensures the existence of a sequence {eK)K converging to such that 

(^f , /„) < fn) + Mf + b /Vf , fn)ds + SK /Vf , l)ds 

Jo Jo 
+ 11 r{x)p{x)a{x)D'^fn{x)pf{dx)ds. (4.2) 

For n > m > 1, splitting M into (— n + m^n — m) and its complement, and using 
Lemma [4. 2 1 we get 

r{x)p{x)a{x)D°' fn{x)iyf (dx) 

<fsupa(x) (snp\D"Mx)\{u^Jn_J + ^^—-{ufA)) . (4.3) 
xm \i,x a(m-l)" J 

Since the sequence (/n)n is non- increasing, {vf , fn) < {vf,fn-m), and there is a 
constant C not depending on n > 2 and m G {2, . . . , n}, and a sequence {r]m)m>2 of 
positive numbers converging to such that 

(i^f , fn) < {y'^. fn) + Mf + C f\uf, fn-n.)ds + {s K + V^) /Vf , 1)^^- (4-4) 

Jo Jo 

Let /if"'" = (sup^<^(z/f- , /„)) and /if- = (sup^<^(z/f^, 1)) which is bounded uni- 
formly in K and t E [0, T] since, using Lemma 14. 1^ 

supE( sup (z/f,l)^) < oo. (4.5) 

-ff te[o,T] 

Assume, first, that < r/ < 1. Observing that, in view of (I2.10p . {M^'"')t < 
CK^~^ 1) + {vf l)'^ds, and using Doob's inequality and (14. Sp . one deduces 

that 

/if'" < /if'" + C /* /if "-'"c^s + + r^^ (4.6) 

for the modified sequences {eK)K and {rjm)m, that still converge to 0, as K and m, 
respectively, grow to oo. Iterating this inequality yields, for j G N*, and m > 1, 

K,n. < ^ x,o-o™ (CT)' OT^-yf/if^. ^ (CT)' 

/ir <2^/^o ~7l~ + + K^K + r]m)2_^—f^ 

l=Q ' \J J- 1=0 ' 



«=L(j+3)/2j ■ yj )■ 
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where we used the monotonicity of /x^'" w.r.t. n, and (14. 5 p to justify the second 
inequahty. 

The random variables {{vq , fn))K converge in law to {C,o,fn), as K ^ oo, and 
are uniformly integrable according to fl3.2p . Therefore, for a fixed n, /x^'" converges 
to (^0) fn), as K oo. Hence, 

limsup/.^'^™< (eo,/L.72j™)e''^+s^P^(«'l)) E 

For m,j large enough, the right-hand-side is arbitrarily small. Since n i— )• /i^'" is 
non- increasing, the proof is complete in the case < rj < 1. 

When 7] = 1, sup„ i?(supj<2.(M^'")t) does not vanish anymore as — j- oo, and 
the proof cannot be concluded in the same way as above. But, taking expecta- 
tions in (^i^ we obtain that holds with supg<j i?((z/f , /„)), E((z/^,/„)) and 
sup^<g E{{h'^ , fn-m)), replacing, respectively, /x^ /Xq and Following the 

previous line of reasoning we obtain that 

lim limsup sup £'((z/f", /„)) = 0. 

n-s>oo K^oo t<T 

Now, in view of f l2.10p . and because < /„, one gets 

(M^'")t <f ! f 2(z/f , /„) + ! Uix + h)MK{x, dh)uf{dx)] ds + Sk, 



with {eK)K tending to as converges to infinity. Since by fl2.6p . for K large 
enough, 

jj„(x + h)M^(:x,dh) = E (/„(i + < /„_i(i) + > 1 

K \ a 



one deduces that limn-s-oo linisup^__^oQ -E((M^'")r) = 0. Taking advantage of (14. 4p 
and Doob's inequality one checks that (14. 6 p holds with Ek replaced by eK,n such 
that lim„_j.oo lim sup^^^^ eK,n = on the right-hand-side. Then one easily adapts 
the end of the proof written for 77 < 1 to obtain the desired conclusion. I 



Remark 4.4 In the case of a mutation kernel with finite second order moment as 
in 0, the true Laplacian being a local operator, (14. 3 p may be replaced by 

The conclusion of Lemma \4-3\ may be derived by a similar argument. 
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Now, let N{dt, dh) be a Poisson random measure on M+ x (—1, 1) with intensity 
c^^ lfep+a • The process Zt = /(ot]x(-i i) h(^N{dt,dh) — jj^tt^^ is a Levy process such 
that, for all p > I, E{\Zi\p) < +00. Let denote the solution of the stochastic 
differential equation 

dX^ = a{X^-)dZu X^ = X (4.7) 

which admits a unique solution when a is Lipschitz continuous according to Theorem 
7, p. 259, in [22]. Let us denote by Pt the associated semigroup defined for all 
measurable and bounded / : M — )■ M by the formula Ptf{x) = E{f{Xf)). The 
infinitesimal generator of the process X^ is 

Lf{x) = j^^^^ [fix + a{x)h) - fix) - f'ix)aix)h) (4.8) 

Lemma 4.5 Assume that a is C"^ with a hounded first-order derivative, and a 
hounded, and locally Lipschitz second-order derivative. For f G C^(M), Ptfix) he- 
longs to C^'^([0,T] X R), and solves the initial-value prohlem 

dtPtfix) = LPJix), it,x) G [0,T] X M 
Pofix) = fix), xeR 



Proof . By Theorem 40, p. 317. in j22], the mapping x Xf is twice continu- 
ously differentiable with first, and second order derivatives solving, respectively, the 
equations 

dd^X^ = a'iX^^)d^X^-dZt, d^X^ = 1, 



and 



dd,,X^ = a'iXnd,,X^-dZt + a"iXnid.Xn'dZt, d^.X^ = 0. 



For g > 1, since EilZ^l'^) < +00, by Theorem 66, p. 346, in |22j, there is a finite 
constant K such that, for any predictable process iHt)t<T, and any t <T, 



E I sup 



s<t 



HrdZr 



< K / Ei\Hs\'')ds. 



(4.9) 



This, combined with the regularity assumptions made on a, and Gronwall's Lemma, 
immediately implies that, for any g > 1, 



sup Ei\d^X^\'^ + n < +CX). 

(i,a;)e[0,T]xIR 



(4.10) 



YoifeCl 



one deduces that the mapping x 1— )■ Ptfix) G C^(M), with d^Ptfix) 



E(/'(Xf)9,Xf) and d^M'ix) = E(r(Xf)(9,Xf)2 + /'(Xf)9,,Xf). Indeed, for 
instance, to differentiate for the second time under the expectation, we observe that. 
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as y tends to x, the random variables {f {X^)dxX^ — f {Xf)dxXf) / {y — x) converge 
a.s. to f" {X^){dxX^Y + f'{Xt)'^xxX^, and by (14.101) are uniformly integrable since 



y-x 



< 



< 



ft VV fl VX 



y-x 



/'(Xf ) - fiXf) 



\\f\ 



y-x 



\dxxX,'\dz+\dxX^ 



y — X 

iinioo 

y~x 



\dxXt\dz. 



Moreover, since /(Xf), /'(Xf)9^Xf , and /"(Xf )(9^Xf)^ + /'(Xf)9^^Xf , are con- 
tinuous w.r.t. X, and right-continuous and quasi left-continuous w.r.t. t, one deduces 
that the mapping (t, x) i— )■ {Ptf{x),dxPtf{x),dxxPtf{x)) is continuous and bounded 
on [0,T] X R. With the upper-bound, 

\PJ{x + a{x)h) - PJ{x) - dxPtf{x)a{x)h\ < \\\dxxPtf\\ood^\x)h\ 

one concludes that the mapping t-)- LPtf{x) is continuous and bounded on 

[0,T] X M. For m > 0, by the Markov property stated in Theorem 32, p. 300, of |22] . 
Pt+ufix) = E{Ptf{XZ)). By Ito's formula. 



P,/(X:) =P,f{x) + / LPJ{X:)ds 



+ 



(P,/(X;_ + a(X,--)/^) - iN{ds, dh) 



dh 



{0,u]x{-l,l) 



Since the last integral is a martingale, it is centered and one obtains that 

Pt+uf{x) - Pi/(X) 



u 



e(^^£ LPj{x:)dsy 



By Lebesgue's Theorem, one deduces that lim„^o+ ^^-f^^^ = LPtf{x). Hence 

{t,x) I—)- Ptf{x) belongs to C^'^([0,T] x M) and solves the initial-value problem 

'dtPtfix) = LPtfix), (t,x) G [0,T] X R 
Po/W = /(x), x e R 



Unfortunately, in the general case needed in the proofs of our theorems, a is 
merely Lipschitz, and Pt-sf{x) is not smooth enough in the spatial variable x. That 

is why, for £ > 0, we set a'^{x) = J^d'{x — y)e~^ dy / \/2'Ke, and define X^'^ as the 
solution of the SDE similar to (14. 7p . but with a replaced by a^. The generator 
of Xl'^ is the operator L'^ defined like L, but with o"^ replacing a. Finally, we set 
P^f{x) = E{f{Xl'^)). Now, we want to know what happens when e tends to 0. The 
next lemma gives the Holder's continuity of dxPt-sf and the order of convergence 
of Pf/(x) toPt/(x). 
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Lemma 4.6 Assume that a is Lipschitz, and f G C|(]R). Then there exists a 
constant C > 0, such that for all t G [0, T], and a;, ?/ G M such that \x ~ y\ < ||(3"||oo, 

\d.PLJ{^) - d.PLJ{y)\ < 

Moreover, 

sup |Pt/(x) - Ptj{x)\ < C^e. (4.11) 

(t,x)e[0,T]xR 

Proof . One has 

+ r E{nxnd.xrd,xr)dz. 



The absolute value of the second term of the r.h.s. is smaller than C\y — x\ since, 
for each g > 1, 

sup sup E(|9^Xf''n < +00, (4.12) 

£>0 {t,x)G[0,T]xR 

because sup^^g 11^^ 'lloo is not greater than the Lipschitz constant of a. To deal with 
the first term, we remark that 

d^^Xr - d.Xr = fa''{X^f){d.X^f - d.Xff)dZs 
Jo 

+ [\a''{x^i')~a''{xyf))d,xyfdz, 

Jo 

Using (14.91) . and the inequality ||(3'^"||oo < one deduces that 

E fsup(9,x:'^ - d^xy'^y] <c f E{{d,xr - d^xy/f)ds 

\s<t J Jo 

^C{x-yl r E ^^Q^x:n\d^Xyn') dzds. 

^ Jo Jy 

In view of (I4.12p . and Gronwall's Lemma, one concludes that 

E (sMd.X:^' - d.Xy/f \ < 9i^^iy)l. (4.13) 

\s<T ) S 

Combining this bound with the inequality 

\E u\xr){d^xr- d.xn) i < ce (i^.x^^ - d.xr\H\d^xr\ + \d.xr\y-^) , 

and Holder's inequality, one easily deduces the first statement of the Lemma. Since 
||<5" — o'^lloo < Cy/^, by a reasoning similar to the one employed to prove fl4.13p . one 
easily checks that E (sup^<'p(X^ — '^)^) < Ce and the second statement follows. 
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5 Proofs of the Theorems 



Proof of Theorem 13.21 The proof of the theorem will be carried out in five steps. 
Let us fix T > 0, and rj < 1. 

Step 1 We first endow Mp with the vague topology. To show the tightness of 
the sequence of laws = C-iy^) in V{p{\Q,T], {Mp,v))), where Mp is endowed 
with the vague convergence topology, it suffices, following Roelly [23], to show that 
for any continuous bounded function / on M the sequence of laws of the processes 
(z/^, /) is tight in D([0, T], M). To this end we use the Aldous criterion [1] and the 
Rebolledo criterion (see [18]) which require us to show that 

supE{ sup ,/)|) < oo, (5.1) 

K te[o,T] 

and the laws, respectively, of the predictable quadratic variation of the martingale 
part, and of the drift part of the semimartingales (z/^, /), are tight. 

Since / is bounded, (15. ip is a consequence of (14. 5p : let us thus consider a pair 
{S, S') of stopping times satisfying a.s. the inequality 0<S<S'<S + 6<T. 
Using (I2.10p and (14. 5p . for some positive real numbers C and C, we get 

E {{M^^f)s' - (M^'Os) < CE (^J^'''' ((z/f , 1) + (z.f , 1)2) ds^ < C'5. 

In a similar way, we show that the expectation of the finite variation part of 
(^^5 /) ~ {^s 1 f) is bounded by C'S. Hence, the sequence = Ciy^) is tight 
inP(D([0,T],(M^,t;))). 

Step 2 Let us now denote by Q the weak limit in ^(©([0, T], (Mi?, y))) of a 
subsequence of (Q^^) which we also denote (Q^)- We remark that by construction, 

sup sup \{ufj)-{uf.J)\<l/K. 

te[o,T] /eL°°(M),||/||oo<i 

Since, for each / in a countable measure-determining set of continuous and com- 
pactly supported functions on M, the mapping z/ supf^pKut, f) — {^t-,f)\ is 
continuous on 3{[0,T], {Mf,v)), one deduces that Q only charges the continuous 
processes from [0, T] into {Mp, v). Let us now endow Mp with the weak convergence 
topology and check that Q only charges the continuous processes from [0,T] into 
{Mp^w), and that the sequence {Q^) in V{p{\Q,T],{Mp^w)) converges weakly to 
Q. For this purpose, we need to control the behavior of the total mass of the mea- 
sures. We will employ the sequence (/„) of smooth functions introduced in Lemma 
14.31 which approximate the functions l{|a;|>n}- For each G N, the continuous and 

compactly supported functions (/„,; '= fn{^ — fi))im increase to /„, as Z — oo. Con- 
tinuity of the mapping z/ i— )■ sup^^p{i't, fn,i) on B){[0,T], {Mp,v)), and its uniform 
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integrability deduced from f l4.5p . imply the bound 

( snp{ut Jr,i)) = lim E (snp{uf J^A < liminf E (snp{uf,U)] . 

\t<T / -R'^oo ) K^co \t<T J 

Taking the hmit, / — )■ oo, in the left-hand-side, in view of the monotone convergence 
theorem and respectively, (14 .Sp and Lemma [4.3^ one concludes that for = 0, 

^« f sup(i/i, 1)^ = E« f sup(i/i, h)] < +00 (5.2) 

\t<T J \t<T J 

and for general n, 

lim E« fsup(i/t,/„)) =0. (5.3) 

\t<T J 

As a consequence one may extract a subsequence of the sequence (supj<'r(t't, /„))„ 
that converges a.s. to under Q, and the set {i't)t<T is tight Q-a.s. Since Q only 
charges the continuous processes from [0,T] into {Mp,v), one deduces that Q also 
only charges the continuous processes from [0, T] into {Mp,w). 

Let X denote a process with law Q. According to Meleard and Roelly [2T], 
to prove that the sequence (Q^) converges weakly to Q in P(D([0,T], {Mf,w)), it 
is sufficient to check that the processes {{i^^ , 1) = ((^'t^, i))t<T)K converge in law 

to {X,l) =^ {{Xt,l))t<T in D([0, T],M). For a Lipschitz continuous and bounded 
function F from D([0, T], M) to M, we have 

lim sup \E{F{{u^, 1)) - F((X, 1))| < lim sup lim sup \E{F{{u^, 1)) - F((z/^, 1 - /„)))| 

+ lim sup lim sup |S(F((z/^, 1 - /„)) - F((X, 1 - /„)))| 
+ lim sup |i?(F((X, 1 - /„)) - F((X, 1)))|. 

n— >oo 

Since |-F((z^, 1 — /„))— F(z/, 1))| < C sup^<'p(t'j, /„), Lemma lOl and (15. 3 p respectively 
imply that the first and the third terms in the r.h.s. are equal to 0. The second term 
is in view of the continuity of the mapping z/ i— (z/, 1 — /„) in D([0, T], [Mp, w)). 

Step 3 Recall that the time T > is fixed, and < 1] < 1. Let us check that, 
almost surely, the process X solves (13. 3p . By (15. 2p . for each T, supjg[o,T] (^t^ 1) is 
finite a.s. Now, we fix a function / G C^(M), compactly supported if a < 1 (in this 
case the extension of (13. 3p to any function / in follows by Lebesgue's theorem), 
and a t < T. 

For u G D([0,T], (Mf,w)), denote 



= {ut, f) - (z/o, /) - / / (6(x, V * Us{x)) - d{x, U * Us{x)))f{x)us{dx)ds, 
Jo Jr 

$2(^) = _ f f a{x)D''f{x)v,{dx)ds. (5.4) 
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We must show that 

EQ{\<fl{X) + ^',iX)\)=0. (5.5) 
By (12.91) . we know that, for each K, 

where 

ft f 



^f^(z/^)=- / p{x){K'^r{x) + b{x,V*uf{x))) 
Jo Jr 

X (^j f{x + h)MK{x,h)dh- f{x)^uf{dx)ds. 
Moreover, (14.51) imphes that for each K, 
E(|Mf'^P) =E((M^'/),) <^i? A) + {uf AY] d.^ < 



K ' 
(5.6) 

which goes to as iiT tends to infinity, since < < 1. Since, by Proposition 13. 
there exists a deterministic sequence (e^), converging to as — )■ oo, such that 



one deduces that 



|v[/?'^(z.^)-Vl/2(^^)|<e^ [\u^<,l)ds, 

Jo 



K 



Since X G C([0,T], {Mf,w)) and / G C^{R), due to the continuity of the param- 
eters, the functions and are a.s. continuous at X. Furthermore, for any 
z/GD([0,T],Mf), 

+ < Cf ({ut, 1) + {uo, 1) + /" 1) + {Us, Ifds] . (5.7) 



Hence, in view of (14. 5p . the sequence (^'^(z/^) + '^^{i'^))k is uniformly integrable, 
and thus 

\imE + = E il^liX) + vl/2(x)|) = 0, (5.8) 

which concludes the proof of the first part of Theorem 13.21 

Step 4 In this step we will prove part (ii) of Theorem 13.21 which asserts unique- 
ness of the solution of (13. 3p under the additional assumption that a is Lipschitz. 
According to (13. 5p . one has 

a(x)D-/(x) = Lfix) + [ (fix + a{x)h) - f{x)) '^^ 

JR\(-1,1) 
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where L has been defined in (14. 8p . It is easy to prove that if ,^ is a solution 
of (13.31) satisfying supjgjo,r] (65 1) < then, for each test function ipt{,x) = 'i/j{t, x) G 
Cl'^{R+ X M), one gets 



= (ecV-o) + / / ids^Pis^x) + aix)D''^ljsix))Udx)ds 
Jo Jr 

+ / {b{x,V * ^s{x)) - d{x,U * is{x)))ip{s,x)^s{dx)ds 
Jo Jr 

= (^0,^0)+ / / {dstlj{s,x) + LilJs{x))^s{dx)ds 



+ / / ({b{x,V*U^))-d{x,U*Ux)))^{s,x) 
Jo Jr V 

+ / ii^is,x + a{x)h) -i;{s,x))—-—]^sidx)ds. (5.9) 

Let t G [0, T] and / G C^(M). We would like to choose ipis, x) = Pt-sf{x), where 
Ptf{x) = -E(/(Xf)), and Xf is the unique solution of the stochastic differential 
equation ( 14. 7p . so that the second term in the right-hand-side above vanishes. In 
view of Lemma 14.51 this is immediately possible if a is assumed to be C^, with a 
bounded first order derivative, and a bounded and locally Lipschitz second order 
derivative. 

Unfortunately, when a is merely Lipschitz continuous, Pt-sf{x) is not smooth 
enough in the spatial variable x. That is why, for e > 0, we set a^{x) = f^^cr^x — 

?/)e~2F^^^j and define X^'^ as the solution to the SDE similar to (14. 7p but with a 
replaced by a^. The generator of X^'^ is the operator defined like L but with 
0"^ replacing a and we set P^f{x) = E{f{Xl'^)). According to Lemma [4.5[ for the 
choice ipi^s^x) = Pj'Lg/(a;), equation (15. 9p takes the form 



(6,/)=(eo,P//)+ r / {L-LnPlJ{x)Ud: 
Jo Jr 



x)ds 



+ / / {h{x, V * Ux)) - d{x, U * Ux)))Plsf{^) 



+ 



/ {PlJ{x + a{x)h) - PlJ{x)) -^)^,(dx)ds. (5.10) 
Jr\(-i,i) I^^I / 



,1) 

One now whishes to make e tend to 0. Since, 

. „ „ .X . dh 

'(-1,1) Ja~{x)h 

Lemma [4.61 and the estimation ||(3" — (3"^||oo < Cv^j imply that 



(L - L^)P!_Jix) = [ r^^^' {d.PlJ{x + y)- d^PlJix)) dy- 

J(-l,l) Ja~(x)h 



|(L - L^)PlJ{x)\ < C '"^^^ / < 



£^4 |/i|«/2 
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Letting e — )■ in (15.101) . one concludes using (14.111) that 



(6,/)=(^o,Pi/) 



+ 



(-1,1) 



{Pt^sf{x + cr{x)h)-Pt-sf{x)) 



dh 



\h\ 



l+a 



^s{dx)ds. (5.11) 



Next, we consider the variation norm defined for /ii,/i2 ^ as follows: 



l/ii - /U2I 



sup 



(5.12) 



By the Jordan-Hahn decomposition of the measure fii—fi2, there exists a Borel subset 
y4 of M such that — /i2|| = (yUi — fJ'2, 1a — 1r\a)- In view of the inner regularity 
of the measure /ii + fi2 there exists a closed set 5 C A such that {fii + fi2, A \ B) is 
arbitrarily small. The function fk{x) = (1 — kd{x, B)) V (—1) is Lipschitz continuous 
and tends to 1b — 1k\b, as k tends to 00. By Lebesgue's theorem, (/ii — /i2, fk{x)) 
tends to (/ii — /i2, 1b — 1r\b), as A; — )■ 00. Since 

1 1 l/Wi - At2| I - (/Ui - /U2, Is - 1r\b) I = 2| (/ii - ;U2, a \ I < 2 + /i2, a \ S) , 

for i? and k well chosen, (yUi — fi2, fk) is arbitrarily close to /i2||- Now, utilizing 
the convolution, fk may be approximated by a sequence of functions globally 
bounded by 1 which converge uniformly on compact sets. Thus one deduces that 

11/^1 - fJ'2\ \ = SUPjgc'2(iR)^ ||/||oo<l I (^1 ~ ^2, /) I- 

Now, we are ready to prove the uniqueness of a solution of (15. lip . For two 
solutions {^t)t>o, and (6)t>o, of (15. lip, such that sup^g^^r] (6 + 6, l) = At < +00, 
and / G C^(M) satisfying the condition ||/||oo < 1, one has 



(6 - L f)\< I I [Udx) - Udx)] ({h{x, V * Ux)) - d{x, U * Ux)))Pt-sf{x) 
Jo Jr \ 

dh 



+ I {Pt.sf\x + a{x)h)-Pt_J{x)) 
i\h\>i 



ds 



+ 



+ 



Udxmx, V * Ux)) - Kx, V * Ux))) 



ds 



Udx){d{x, U * Ux)) - dix, U * Ux)))Pt-sfix) 



ds. 



(5.13) 



Since 



,0 < 1, then ||-Pt_s/||oo < 1 and, for all x G M, 

\{b{x, V * Ux)) - d{x, U * Ux)))Pt-sf{x)\ <b + d{l + UAt), 
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and 



/ (P,_ J(x + a{x)h) - Pt-sfix))-^ 
J\h\>i 1^1 + 



< 



a 



Moreover, b and d are Lipschitz continuous in their second variable with respective 
constants Kb, and Kd. Thus we obtain from f l5.13p that 



(6-6,/>l < 



b + d(l + UAt) + - + KbArV + K^AtU 
a 



\\is-U\ds. (5.14) 



Taking the supremum over Cl functions / bounded by 1, one obtains 



116-611 < 



h + d{l + UAt) + - + K^AtV + K^AtU 
a 



116 -61 Ms, 



and uniqueness follows by an application of Gronwall's Lemma. 



Step 5 In the final step we shall prove the existence of a density claimed in 
Part (in). The image of the Poisson random measure N on M_|_ x [—1,1] with 
intensity dt-p^pp^, by the mapping {t,h) (t, ^p^), is a Poisson random measure 

on M+ X R \ [— ^] with intensity dtdz since, for z = one has dz = jj^tt^- 

Let us denote by p,{dt, dz) the associated compensated measure. The stochastic 
differential equation (14. 7p can now be written in the form. 



(0,t]xM\[-l 1] 



c(X^-, z)fi{ds, dz), 



for c{x,z) = a(x) x . When the strictly positive function a E C|, i.e., it is 

(a\z\)a 

bounded together with its derivatives up to order 3, one may apply Theorem 2.14, 
p. 11, to deduce that, for t G (0,T], Xf admits a density pt{x,y) with respect to 
the Lebesgue measure on the real line. With flS.lip . one deduces that, for t > 0, C,t 
admits a density equal to 

Pt{x,y)^o{dx) + / / {b{x,V *^s{x)) - d{x,U * ^s{x)))pt-s{x,y)^s{dx)ds 
J«. Jo Jr 

dh 




+ 111 {Pt-s{x + a{x)h,y)-pt-s{x,y))——is{dx)ds. 

'0 JRJ\h\>l \l>'\ 

This completes the proof of Theorem 13. 2[ 



□ 



Let us now turn to the proof of Theorem 13.31 In this case the allometric exponent 

7] = 1. 
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Proof of Theorem 13.31 We will use a method similar to the one employed in the 
proof of Theorem I3.2[ Actually, Steps 1, 2 and 3 are completely analogous and we 
omit them. Thus, we only have to prove the uniqueness (in law) of the solution of 
the martingale problem fl3.7l) - fl3.9p . and the fact that any accumulation point of the 
sequence of laws of is a solution of fl3.7p - fl3.9p . 

Uniqueness. The uniqueness in the general case can be deduced from the special 
case when 6 = ci = by using the Dawson-Girsanov transform for measure-valued 
processes (cf. Evans and Perkins [13] (Theorem 2.3)). Indeed, 



which allows us to use this transform. 

In the case h = d = Q the proof of uniqueness can be adapted from Fitzsimmons 
[H], Corollary 2.23. This proof is based on the identification of the Laplace trans- 
form of the process, using the extension of the martingale problem (13. 8p to functions 
'?/'(s,x) = Pt-sf{x) with bounded functions / (Fitzsimmons [H] Proposition 2.13). 
For iIj{s,x) being in time, and in x, (13. 8p extends with an additional term 
dgip appearing in the drift part. As we have seen above, the function Pt_sf{x) is 
not smooth enough. Thus we firstly apply (13. 8p to ■?/'(s,x) = P^_gf{x), as already 
done in the proof of Theorem 13. 2[ thanks to Lemma 14.51 and we make e tend to 
as in Lemma [4.61 

Identification of the limit. Fnally, let us identify the limit of the sequence 
of laws of as a solution of ([X7D-(IS3D- Write = C{v^), and denote by 
Q a limiting value in P(D([0, T], (M^, w)) of a subsequence (denoted also Q^), 
and by X = {Xt)t>o a process with law Q. Because of Step 4, X belongs a.s. to 
C([0, T], (Mi?, )). We have to show that X satisfies the conditions (13. 7p . (13. Sp 
and (13. 9p . First, note that (13.70 is straightforward from (13.60 and Lemma I4T1 

Next we show that for any function / in C^(]R), compactly supported when 
a < 1, the process M/ defined by is a martingale (the extension to every 

function in Cl is not hard). Indeed, consider < si < ... < s„ < s < t, and 
continuous bounded maps 01, on Mp. Our goal is to prove that, if the function 
\E' from ©([0, T], Mp) into R is defined by the expression 



E 




n 



k=l 




then 



E(^(X)) = 0. 



(5.16) 
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It follows from that 

= ( n MO {Mf - Mf'f] ) = E inO) - Ak, (5.17) 



where 



AK = E(f[MO f I {p{x)b{x,V*u^{x))\ [ {f{x + h)-f{x))MK{x,dh) 
+ p{x)r{x)K [ {f{x + h)-f{x))MK{x,dh)-d{x)D''f{x)]u^{dx)du 



In view of Proposition I3.H Ak tends to zero, as K grows to infinity. Applying 
Lemma [4.H for p = 3, we see that the sequence (|\E'(z/^)|)i^ is uniformly integrable, 
so that 

limE{\^iu^)\)=EQi\^{X)\), (5.18) 

since the function ip is continuous a.s. at X. Collecting the previous results allows 
us to conclude that fl5.16p holds true, and thus M'^ is a martingale. 
Finally, we have to show that the bracket of M-^ is of the form 



(M^)i = 2 / r{x)f{x)Xs{dx)ds. 
Jo Jr 

To this end, we first check that 

Ni = {Xt, fY - (Xo, JY- [ f 2r{x)f{x)X,{dx)ds 

Jo Jr 

-2 [ (XsJ) [ fix)[bix,V*Xsix))-d{x,U*Xsix))]Xsidx)ds 



-2 [ (X,,/) / a{x)D^f{x)X,{dx)ds (5.19) 
Jo Jr 

is a martingale. This can be done exactly as for M/, using the semimartingale 
decomposition of {uf, f)^, given by (12.81) with 0(z/) = (z/, /)^, and applying Lemma 
14.11 with p = 4. On the other hand, Ito's formula implies that 



(Xi,/)2-(Xo,/)2-(M^)i-2 / (X,,/) / a{x)Dy{x)X,{dx)ds 

Jo Jr 

-2 [ (X,,/) / f{x)[bix,V * Xsix)) - d{x,U * Xsix))]Xsidx)ds 
Jo Jr 

is a martingale. Comparing this formula with fl5.19p . we obtain (13.91) . □ 
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6 Concluding remarks 

We have developed models for population dynamics in the context of evolutionary 
ecology permitting heavy tailed distribution of mutations. Depending on the value 
of the allometric exponent t], the continuum (macro) limits of the individual (micro) 
dynamics turned out to be described by deterministic solutions of fractional nonlocal 
reaction-diffusion equations driven by fractional Laplacians (the case < rj < 1), 
or measure-valued nonlinear stochastic superprocesses driven by Levy- stable pro- 
cesses. These limiting models can now be used as approximate objects for numerical 
simulation of evolutionary Darwinian dynamics in presence of non-negligible large 
mutations. Of course, estimators of the relevant parameters of the phenomena under 
study have to be obtained first. 

It the future it may also be worthwhile, from the perspective of practical applica- 
tions, to elucidate the situation where the mutations have distributions intermediate 
between the heavy tailed distributions studied in this paper and the Gaussian dis- 
tributions considered in p]. Such distributions, which can display a multiscaling 
behavior, a-stable type for small mutations, and exponential for large mutations, 
have been recently suggested in the physical and economics literature and studied 
under different names such as truncated Levy, and tempered Levy distributions, see, 
e.g., 0, ED, ES], and [ID]. 
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